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Abstract: Software which corrects the dynamic error of force transducers in impact force 
measurements using their own output signal has been developed. The software corrects the 
output waveform of the transducers using the output waveform itself, estimates its 
uncertainty and displays the results. In the experiment, the dynamic error of three 
transducers of the same model are evaluated using the Levitation Mass Method (LMM), in 
which the impact forces applied to the transducers are accurately determined as the inertial 
force of the moving part of the aerostatic linear bearing. The parameters for correcting the 
dynamic error are determined from the results of one set of impact measurements of one 
transducer. Then, the validity of the obtained parameters is evaluated using the results 
of the other sets of measurements of all the three transducers. The uncertainties in the 
uncorrected force and those in the corrected force are also estimated. If manufacturers 
determine the correction parameters for each model using the proposed method, and 
provide the software with the parameters corresponding to each model, then users can 
obtain the waveform corrected against dynamic error and its uncertainty. The present status 
and the future prospects of the developed software are discussed in this paper. 
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1. Introduction 

Force transducers are widely used in both the industrial and research areas. However, at present 
only static force calibration methods using static standard forces are available. As a result of this fact it is 
very difficult for users to estimate the dynamic errors of force transducers and correct those errors. 

A method, i.e., the Levitation Mass Method (LMM), has been proposed for dynamic force 
calibration [1,2]. In the LMM, varying forces, such as an impact force, an oscillation force and a step 
force, are applied to the transducer as the inertial forces of a mass supported using an aerostatic linear 
bearing, i.e., the moving part of the bearing. Then, the forces measured as the inertial force are 
compared with the output signal of the transducer. In the method, only the Doppler frequency shift 
of the laser light reflected on the mass is accurately measured using an optical interferometer. Then all 
the other quantities, such as velocity, acceleration, displacement and inertial force, are calculated from 
the frequency. 

Xu et al. proposed a method in which the inertial mass of a part of the transducer itself is supposed 
to be the error source and its value is estimated from the free oscillation of the transducer after 
applying an impact load [3]. However, in the method, the real force during the impact is unknown and 
that is not used in the calibration. 

Using the LMM, Fujii et al. experimentally showed that the error in impact force measurement is 
almost proportional to the acceleration at the sensing point of an S-shaped strain-gage force transducer, 
and further, it can be explained as the effect of the inertial force of a part of the transducer itself. The 
difference between the static and dynamic characteristics of the force transducer is well explained 
as the inertial force of a part of the transducer itself [4]. This result suggests that the dynamic 
measurement errors can be corrected by measuring the acceleration at the sensing point. 

Based on this result, a force sensor, which consists of a strain-gage force transducer and an 
accelerometer, is proposed [5]. The effect of the inertial mass of a part of the sensor is corrected using 
the acceleration measured using the accelerometer. 

Then, it is shown that the dynamic errors of force transducers can be corrected using their own 
output signal [6]. It is experimentally confirmed that the acceleration at the sensing point of the 
transducer is almost proportional to the differential coefficient of the second order of the output signal 
of the transducer itself. This can be understood by the fact that the deformation of the transducer 
is almost proportional to the force applied. However, only one transducer is used in the experiment. 

In this paper, the validity of the proposed method for correcting the dynamic error and for 
evaluating its uncertainty is experimentally evaluated using three force transducers of the same model. 
Based on the experimental results, software, which corrects the dynamic error of force transducers 
using their own output signal, has been developed. We suggest that if the sensor manufacturers provide 
the software with sets of the parameters corresponding to their models, then the users can obtain the 
corrected waveforms and reduce the uncertainty in dynamic force measurements. 
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2. Experiment 

Figure 1 shows schematic diagram of the experimental setup for measuring the impact response of 
force transducers. The transducer under test is fixed on a large metal base of 10 kg. An aerostatic 
linear bearing is used in order to obtain a linear motion with sufficiently small friction acting on the 
mass, i.e., the moving part of the bearing. The initial velocity is given to the moving part by hand. 

Figure 1. Experimental setup. Code: LD = Laser diode, PD = Photo detector, CC = Corner 
cube prism, PBS = Polarizing beam splitter, NPBS = Non-polarizing beam splitter, 
GTP = Glan-Thompson prism, ADC = Analog-to-digital converter, DAC = Digital-to-analog 
converter, DSR = Dynamic strain recorder. 
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A rubber block of sufficiently small mass is attached at the sensing point in order to adjust 
the steepness of the impact. The velocity of the mass vi and the velocity of the sensing point of 
the transducer V2 are measured using two optical interferometers named Interferometer- 1 and 
Interferometer-2, respectively. The force acting on the sensing point is obtained by the equation of 
motion, F = ma, where m is the mass and a is the acceleration. 

The mass of the moving part including the cube corner prism (CC) and the extension rod, M\, 
is approximately 2.653 kg. The mass of the contact point of the transducer, including the CC and their 
base plate, M 2 , is approximately 0.082 kg. 

Three conventional S-shaped strain-gage type transducers (model: DB-200N, capacity: 200N 
manufacturer: Showa Measuring Instruments Co. Ltd., Tokyo, Japan) — Transducer- 1, Thransducer-2 
and Transducer-3 — are used in the experiment. The year of manufacture of Transducer- 1 is 2000, and 
that of the other two transducers is 2008. 

A Zeeman-type two-frequency He-Ne laser is used as the light source of the optical interferometers. 
The interferometers have three photo-detectors: PD0, PD1 and PD2. The frequency difference between 
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the two orthogonal polarization states emitted from the laser, / rest , is monitored using a Glan-Thompson 
prism (GTP) and the first photo -detector, PDO. 

The total force acting on the moving part of the aerostatic linear bearing, F mass , is calculated as the 
product of its mass, M\, and its acceleration, a\, as follows: 

F ma ss = Miai (1) 

In the measurement, the total force acting on the mass, Fmass, is considered to be the same as the 
force acting on the mass from the force transducer being tested, since the frictional force acting on the 
mass is negligible [7]. The acceleration is calculated from the velocity of the levitated mass. The 
velocity of the mass, i.e., of the moving part of the aerostatic linear bearing, vi, is measured as the 
Doppler shift frequency, /boppieri, which can be expressed as follows: 

Vl = ^air(/bopplerl)/2 (2) 
/bopplerl = — (/beatl — /rest) (3) 

where ^ a ir is the wavelength of the signal beam under the experimental conditions, and/ b eati is the beat 
frequency, which is the frequency difference between the signal beam and the reference beam and 
appears as the beat frequency at PD1. The position of the mass x\, the acceleration of the mass a\ and 
the force acting on the mass, F mass , are numerically calculated from the velocity. 

On the other hand, the velocity of the sensing point of the force transducer, V2, is measured as the 
Doppler shift frequency, / D o PP ier2, which can be expressed as follows: 

V2 = ^air(/boppler2)/2 (4) 
/boppler2 = — (/beat2 — /rest) (5) 

The beat frequency /b ea t2 is measured using PD2. The position X2 and the acceleration a 2 of the 
actuator are numerically calculated from the velocity v 2 . 

For simultaneous measurement, three frequency counters and the dynamic strain recorder (DSR) are 
triggered by a single signal originated by a light switch composed of a laser diode (LD) and a PD. 

The force measured by the transducer, F tr an S , is calculated using the output signal of the transducer 
Virans stored by a dynamic strain recorder (DSR; model: DC-204R, Tokyo Sokki Kenkyujo, Tokyo, 
Japan) and its static calibration result. The force calculated using the output signal of the transducer 
and its static calibration result, F ttms is compared with the force measured as the inertial force of the 

maSS, Fmass- 

In the experiments, one set of experiment consists of 20 impact measurements. The first set of 
the experiments using Transducer- 1, Set_Al, is conducted for evaluating the parameters for dynamic 
error correction. Then, additional three sets of experiments using all the three transducers are 
conducted for evaluating the validity of the estimated parameters. These additional three sets of the 
experiments using Transducer- 1, Transducer-2 and Transducer-3, are described as Set_Bl, Set_B2 and 
Set_B3, respectively. 
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3. Estimation of the Parameters for Dynamic Error Correction 

Figure 2 shows the result of a single set of the impact tests of Transducer- 1 . Figure 2 shows 
the force measured by the force transducer Ftrans, and the force acting on the mass F mass and then- 
difference Fdiff: 

F&fif = Frans ~~ ^mass (6) 

Figure 2. The force measured with the transducer Ftrans, the force measured with the 
interferometer F mass and the difference between the forces. 
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The figure shows the electric response of the force transducer under the impact load. The difference 
between F tmns and F mass is derived mainly from the difference between the static and dynamic 
characteristics of the transducer. The root mean square (RMS) value of the difference between F tmDS 
and F mass during the collision period, 0 ms < t < 14 ms, was approximately 6.2 N. Figure 3 shows the 
relationship between the second time derivative of the force measured using the force transducer, 
d 2 F tmn Jdt 2 , and Fdiff, for Set_Al. 

Figure 3. Relationship between the differential coefficient of the second order of the force 
measured with the transducer and the difference between F ttans and F ma ^. 
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The plots show a linear relationship between d F Um Jdt and Fdiff. The solid line in Figure 3 shows a 
regression line of the following equation. 

F reg = C(d 2 F trans /d? 2 ) (7) 
The parameter for dynamic error estimation C is obtained as the slope of the regression line. Here, 

~~ 7 2 

the regression line is determined by a least-square method and resulted in C = -2.49 x 10 s . The 
measured plots coincide well with the regression line. 

Figure 4 shows the difference between F^ns and F mass and inertial forces of the transducer estimated 
with the parameter C, F reg . The measurement error, Fdiff, is almost the same as the result of the 
regression analysis, F reg . From this result, the corrected force, F cor rected, can be calculated using 
following equation: 

^corrected = ^trans ~~ ^reg = ^trans ~~ C (d^trans/d? 2 ) = Ftrans _ (— 2.49 x 10 ? ) {^F tvim Jdt 2 ) (8) 

Figure 4. Difference between F trans and F mass and estimated inertial force of the transducer. 
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Figure 5 shows RMS values of the dynamic errors Fditf and that of residual errors after 
dynamic correction iv s : 

^res = -^corrected — F mass (9) 

for Set_Al. The RMS values of F&ff and F res seem to be almost square functions of the maximum 

force -f\rans,max- 

Solid and dashed curves in Figure 5 are regression curves of the RMS values of Fam and F res , 
respectively. The RMS value of the differences between Fdiff and its regression line and that between 
F res and its regression line are 1.34 x 10 N and 3.58 x 10 N, respectively. This results show that the 
RMS values of the F dif f and F res can be estimated from the F tranSjmax . 
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Figure 5. RMS valuses of the dynamic measurment error Fdiff = Ftrans 
values of the residual error of the corrected values F res = Fcon-ected _ ^mass- 
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4. Evaluation of the Validity of the Estimated Parameters 

To evaluate the validity of the dynamic error correction using the parameter C estimated using 
Transducer- 1, the dynamic error correction is applied to the other sets of experiments performed by the 
three transducers, i.e., Set_Bl, Set_B2 and Set_B3. Figure 6 shows relations between the full width 
at half maximum (FWHM) of F ttms and the maximum force, F tranSjmax for all the four sets of the 
experiments. All the sets show the similar relation. 

Figure 6. Pulse widths of F tmn% against various maximum forces measured by using the 
three transducers. 
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Figure 7 shows the RMS values of the dynamic error F^m and that of residual error after dynamic 
correction F ms , about Set_Bl, Set_B2 and Set_B3. The RMS values of F res are effectively smaller than 
that of Fdiff. The correction with the parameter C has enough of a result even though Transducer- 1 is 
8 years older than the other two transducers. 
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Figure 7. RMS values of F di f f = Ftrans ~ ^mass and F ies 
20 sets of impact test with the three transducers. 
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5. Estimation of Uncertainty 

The standard uncertainty w trans in the un-corrected force F trans is estimated as follows: 

/ 2 " 2 
"trans = V "diff + M LMM 



(10) 



where Wdiff is the RMS value of Fdiff during the impact shown in Figure 5. In the experiments, wlmm is 
estimated to be 1.2 x 10~ 4 N [7] and is negligible. As shown in Figure 5, the w dif f can be assumed as 
a square function of the peak value of the force F tt ans,max as follows: 

'2 



"trans ~ ^S,l ' ^trans,max + ^S,2 ' ^transjnax + ^S,3 



(ID 



where the constants Es,i, £s,2 and £s,3 are parameters for estimating the dynamic error in the un-corrected 
force Ftrans calculated based on the conventional static calibration. 

As same as the Wtrans, the standard uncertainty w CO n-ected in the corrected force F corrected is estimated 
as follows: 



u 



corrected 



2 2 
res + "LMM 



(12) 



where w res is the RMS value of the residual error F res during the impact shown in Figure 5. The w res can 
be assumed as a square function of the peak value of the force F tranSjmax as follows: 



"corrected - ^D,l " ^trans,max + ^D,2 " ^trans,max + ^D,3 



(13) 



where the constants £b,i, £b,2 and Eu,3 are parameters for estimating the residual error in the corrected 
force F corre cted estimated by means of the proposed dynamic correction method. 
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6. Correction Software 

From the knowledge obtained from the above experiments, software that corrects the dynamic error 
in the impact force measurement is developed [8]. It works in the Microsoft Windows environment. 
Figure 8 shows a screenshot of the software. 



Figure 8. Screenshot of the correction software. 
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This software processes the input waveform measured by the transducer based on the static 
calibration result, F trans , into (a) waveform of dynamic error F mg [= C(d 2 F trans /d? 2 )]; (b) waveform of 
corrected force F corrected (= F trans - F reg ); In addition, this software calculates: (c) the standard 
uncertainty w t rans in the un-corrected force F^ans! and (d) the standard uncertainty ^corrected in the corrected 

force F correc ted- 

The software outputs them to a output file and shows them to the monitor window. The format of 
the input waveform file is a text format described the sampling interval in time and the sequence 
of force in newton calculated based on the static calibration result. The format of the output data file is 
a CSV format. The input data and the output data of the developed software are as follows, 
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(A) Input data, which should be manually input from the operational window. 

(1) Name of input text- file, in which the data set of the force calculated based on the static 
calibration results F UsSii {i) with the measurement time t(i) contains: [t(i), F trans (i)] (0 ^ i < N - 1). 

(2) Parameter for dynamic error correction: C. 

(3) Parameters for estimating the dynamic error in the un-corrected force F Xvan;i calculated based on 
the conventional static calibration: Es,\, £s,2 and £s> 

(4) Parameters for estimating the residual error in the corrected force ^corrected estimated by means 
of the proposed dynamic correction method: 2s D ,i, £b,2 and £b> 

(B) Output data, which is written in a CSV-file and shown on the display: 

(1) Waveform of F ttans , F ieg [= C(d 2 F tm Jdt 2 )] and F corrected (= ^trans - ^ reg )- 

(2) Parameters for F trans : the maximum value i 7 t rans,max, the pulse width FWHM tra ns and estimated 
error Mtrans- 

(3) Parameters for F cor rected: the maximum value Fcon-ected.max, the pulse width FWHM corT ected and 
estimated error ^corrected- 

(4) RMS value of the difference between F trans and ^corrected- 

7. Conclusion 

The validity of the proposed method for correcting the dynamic error in the impact force 
measurement is experimentally evaluated using three force transducers of the same model. Based on 
the experimental results, software, which corrects the dynamic error of force transducers using their 
own output signal, has been developed. Because there are no dynamic calibration methods for force 
transducers available at this moment, we believe that the work described in this paper is very valuable. 

There are, however, still residual errors in the corrected force. For example, there is a small error 
factor which is assumed to result from the viscosity of the internal structure of the transducer. For 
more accurate correction and more accurate estimation of the uncertainty of measurement, our group is 
working on more accurate measurement of the force with the LMM and investigating the causes of the 
residual errors. Additionally, we will apply this correction method to transducers of other types and 
other models with wider ranges of forces in the future. 

In this paper, the dynamic error correction is performed considering only with the inertial force of a 
part of the transducer itself. If a certain part is attached to the sensing point of the transducer and the 
force acting to the opposite side of the attached part should be measured, then the correction parameter 
should be modified for the mass of the attached part. 

If the sensor manufacturers were to determine the correction parameters for each model using the 
proposed method; and provide the software with the parameters corresponding to each model; then the 
users could obtain waveforms corrected against dynamic errors and their uncertainty. We believe that 
in the near future; all the manufacturers of force transducers will provide similar software and the 
parameters corresponding to their products, so that transducer users will be able to perform dynamic 
error corrections and remove their uncertainty easily. 
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